Chaotic front dynamics in semiconductor superlattices 
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We analyze the dynamical evolution of the current and the charge density in a superlattice for 
fixed external dc voltage in the regime of self-sustained current oscillations, using a microscopic 
sequential tunneling model. Fronts of accumulation and depletion layers which are generated at 
the emitter may collide and annihilate, thereby leading to a variety of different scenarios. We find 
complex chaotic regimes at high voltages and low contact conductivities. 
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I. INTRODUCTION 

Electronic transport in semiconductor superlattices 
(SL) is known to show strongly nonlinear spatio-temporal 
dynamics. Either self-sustained current oscillations 
jl], ||, ||, [|, H or a sawtooth-like current-voltage charac- 
teristic with many branches associated with static field 
domains ||, [t| ||, [j| have been found. For a recent review 
sec [|l0|. Under time-dependent external voltage condi- 
tions, superlattices exhibit a rich menagerie of complex 
behavior including ac driven chaos and switching sce- 
narios between multistable states. This was studied re- 
cently experimentally [jn], |l2|, [H^, [IJ, [l5| and theoreti- 
cally 

In this paper we present simulations of dynamic sce- 
narios for superlattices under fixed time-independent ex- 
ternal voltage in the regime where self-sustained dipole 
waves are spontaneously generated at the emitter. The 
dipole waves are associated with traveling field domains, 
and consist of electron accumulation and depletion fronts 
which in general travel at different velocities and may 
merge and annihilate. We find that depending on the 
applied voltage and the contact conductivity, this gives 
rise to various oscillation modes and self-synchronization 
effects as well as different routes to chaotic behavior. 

A similar scenario of merging fronts was recently found 
in the context of a spatially continuous model describing 
bulk impurity impact ionization breakdown Jig] . It is also 
reminiscent of patterns of temperature pulses in globally 
coupled heterogeneous catalytic systems, e.g. p9| . 



J m ^ m+ i(F m , n m , n m+ i) from well to to well m + 1 de- 
pends only on the electric field F m between both wells 
and the electron densities n m and n m +i in the wells (in 
units of cm -2 ). For details of the microscopic calculation 
of J ra _» m+ i we refer to the literature [lTj]. A typical 
result for the current density vs electric field character- 
istic is depicted in Fig. [j] in the spatially homogeneous 
case, i.e. Ti m = Ti m -L\ = Nd, with donor density ATrj. 
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FIG. 1: Current density vs electric field characteristic at the 
emitter barrier (straight line) and between two neutral wells. 
The Ohmic conductivity of the emitter is a = 0.5 fU 1 ™," 1 . 
The inset shows the front velocity vs. current density for elec- 
tron depletion and accumulation fronts. Jo and Jt denote 
the current for dipole and tripole propagation, respectively. 



II. THE MODEL 

Weakly coupled superlattices are successfully described 
by a one dimensional sequential tunneling model for elec- 
trons [2C], |TJ. In the framework of this model elec- 
trons are assumed to be localized at one particular well 
and only weakly coupled to the neighboring wells. The 
tunneling rate to the next well is lower than the typ- 
ical relaxation rate between the different energy levels 
within one well. The electrons within one well are then 
in quasi-equilibrium and transport through the barrier 
is incoherent. The resulting tunneling current density 



In the following we will adopt the total number of elec- 
trons in one well as the dynamic variables of the system. 
The dynamic equations are then given by the continuity 
equation 

e ~~TT~ = Jm-l^m - for TO = 1, ... iV (1) 

at 

where N is the number of wells in the superlattice. 

The electron densities and the electric fields are cou- 
pled by the following discrete version of Gauss's law 

£r^o{F m - F m -i) = e(n m - N D ) for to = 1, ... N, (2) 
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where e r and eo are the relative and absolute permittivi- 
ties, e < is the electron charge, and F and F N are the 
fields at the emitter and collector barrier, respectively. 

The applied voltage between emitter and collector 
gives rise to a global constraint 



U 



N 



(3) 



where d is the superlattice period. 

The current densities at the contacts are chosen such 
that dipole waves are generated at the emitter. For this 
purpose it is sufficient to choose Ohmic boundary condi- 
tions: 



Jo- 



aF 



n N 



Jn^n+i — ctFn—— 



(4) 
(5) 



where a is the Ohmic contact conductivity, and the factor 
njsf/Njj is introduced in order to avoid negative electron 
densities at the collector. Here we make the physical 
assumption that the current from the last well to the 
collector is proportional to the electron density in the last 
well. It is in principle possible to calculate the boundary 
conditions using microscopic considerations |2^, |2^] , but 
the qualitative behavior is not changed. 



III. NUMERICAL RESULTS 

In our computer simulations we use an N = 100 super- 
lattice with Alo.3Gao.7As barriers of width b = 5nm and 
GaAs quantum wells of width w = 8nm, doping density 
No = 1.0 x 10 11 cm -2 and scattering induced broadening 
r = 8meV at T = 20K. The contact conductivity a is 
choosen such that the intersection point with the homo- 
geneous current density vs field characteristic in Fig. [j] 
is at a current value at which no stationary field domain 
boundaries exists. By this configuration, accumulation 
and depletion fronts are generated at the emitter. For 
large values of a, e.g. a — 1.3 f2 -1 m -1 , we find that 
those fronts form a dipole, i.e. a traveling field domain, 
with leading electron depletion front and trailing accu- 
mulation front (cf Fig. || (a)). The dipole traverses the 
sample at almost constant velocity and constant current 
which shows up as a plateau in the time trace of the cur- 
rent. The current Jd and the front velocity in this dipole 
propagation mode are given by the intersection point of 
the velocities of accumulation and depletion front (cf. in- 
set of Fig. |l|). As the leading depletion front reaches the 
collector, there is no speed constraint on the remaining 
accumulation front. As it accelerates, the current rises, 
and a new dipole is generated at the emitter. Now for a 
short time there are two accumulation and one depletion 
fronts in the sample. In order to fulfill eq.(^) for fixed U 
the depletion front has to assume twice the velocity of the 
accumulation fronts. This constraint fixes the current to 



Jt (cf. inset of Fig. |l|) during this tripole regime. Note 
that in the current time trace the fast small-amplitude 
oscillations (due to well-to- well hopping of depletion and 
accumulation fronts in our discrete model) in the dipole 
and tripole regime are not resolved temporally. 





150 160 170 180 190 200 

Time [ns] 

FIG. 2: Dynamic evolution of charge density and current 
I for contact conductivities (a) a = 1.3f2 _1 m _ , (b) a — 
O^fi-V" 1 and (c) a = 0.55fi~ f m _1 for fixed bias U = 1.0V. 
Light and dark regions denote electron accumulation and de- 
pletion fronts in the space-time plots of the charge densities, 
respectively. 



At lower contact conductivity a — 0.6 fi m , we 
find that instead of a dipole only a depletion front is gen- 
erated at the emitter, as the old depletion front reaches 
the collector (cf. Fig. || (b)). Now for a short time a 
dipole with leading accumulation front exists. The cur- 
rent is fixed again by the constraint of equal velocities of 
accumulation and depletion front, as for the dipole with 
leading depletion front. At the time the old accumula- 
tion front reaches the collector a new accumulation front 
is generated at the emitter. This process is accompanied 
by a dip in the current time trace. 

For even lower a = 0.55 f2 _1 m _1 the fronts at the 
emitter are generated as dipoles with leading accumula- 
tion and trailing depletion fronts (cf. Fig. || (c)). The 
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FIG. 3: (a) Positions where accumulation and depletion fronts 
annihilate vs voltage at a = 0.5 fi m . The grayscale in- 
dicates high (black) and low (white) numbers of annihilations 
at a given well, (b) Time differences between consecutive 
maxima of the electron density in well no. 20 vs voltage at 
a = 0.5 frV -1 . Time series of length 600ns have been 
used for each value of the voltage. 



velocity of the fronts is again determined by the cur- 
rent and the number of fronts in the sample. Since the 
two types of fronts in general move at different veloci- 
ties, merging and annihilation of accumulation and de- 
pletion fronts may occur, which may lead to complicated 
behaviour including chaos as shown in Fig. ^| (c). 

In order to study the bifurcation scenario leading to 
chaos we now fix the boundary conductivity to a = 
0.5 Q m . Experimentally, low contact conductiv- 
ities may be realized by changing the doping and barrier 
structure of the emitter. The temperature dependence of 
the emitter current may also be exploited. In Fig. || (a) 
a density plot of the positions (well numbers) at which 
two fronts annihilate is shown as a function of the volt- 
age. We see that for low voltage the annihilation takes 
place at one definite position in the superlattice with a 
variation of only a few wells. This distribution broadens 
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FIG. 4: Dynamic evolution of charge density for various volt- 
ages at a — 0.5 f2 _1 m _1 . Light and dark regions denote 
electron accumulation and depletion layer, respectively. 



for increasing voltage in characteristic bifurcation sce- 
narios reminiscent of period doubling, leading to chaotic 
regimes. We note that in the chaotic region periodic win- 
dows exist. Since we are dealing with a discrete system, 
the position of the merging of two fronts is also discrete. 
For more accurate analyses it is convenient to use a con- 
tinuous variable such as the time difference between two 
maxima in the electron density in a specific well. The cor- 
responding bifurcation diagram for well no. 20 is shown 
in Fig. H (b). It exhibits a complex structure. 

The transition from periodic to chaotic oscillations is 
enlightened by considering the space-time plots for the 
evolution of the electron densities at different voltages U 
(Fig. |]). While for U — 0.5V the oscillations are regular, 
we see that at U — 0.7V long and short front patterns 
alternate, which is characteristic for a period doubling bi- 
furcation. Every second pair of fronts travels farther into 
the superlattice before they annihilate. At U — 0.79V 
a further period doubling occurs, as can be most clearly 
seen at the alternating length of the longer front patterns. 
This yields a period-four cycle. Further increase of the 
voltage finally leads to well developed chaos, with front 
patterns of different lengths. For U > 1.8V we find that 
fronts may even traverse the sample, while the dynamics 
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remains chaotic. 



time-independent external voltage conditions, whereas 
previous studies have focussed on ac-driven chaos. 



IV. SUMMARY 

We have shown that complex chaotic spatio-temporal 
scenarios may arise in weakly coupled superlattices under 
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